Adittional levels between Landau bands due to vacancies in graphene: 

towards a defect engineering 
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We describe the effects of vacancies on the electronic properties of a graphene sheet in the presence 
of a perpendicular magnetic field: from a single defect to an organized vacancy lattice. An isolated 
vacancy is the minimal possible inner edge, showing an ant idot like behavior, which results in an extra 
level between consecutive Landau levels. Two close vacancies may couple to each other, forming a 
vacancy molecule tuned by the magnetic field. We show that a vacancy lattice introduce an extra 
band in between Landau levels with localization properties that could lead to extra Hall resistance 
plateaus. 

PACS numbers: 73.43.-f, 73.23.-b, 73.63.-b 



I. INTRODUCTION 

Graphene is one of the topical stages for research in 
condensed matter in the last few years. Since the semi- 
nal work describing the obtention of single graphene lay- 
ers [1^ , a great amount of surprising results was obtained 
121 [3] . Of special interest for the present work is the obser- 
vation of the quantum Hall effect in graphene [4, 5 and 
the related discussion in the origin of the disorder effects 
and the role of vacancies in this context. Vacancies are 
natural defects in the graphene lattice and can also be 
externally induced by ion-beam irradiation. Among the 
theoretical studies dealing with the effects of vacancies 
on the electronic structure of graphene, the main focus 
until now is on situations in the absence of magnetic field 
[HI |7l El E] • Within the experimental framework, images 
identifying atomic vacancies in graphene have been re- 
ported [10 . However, although the subjacent motivation 
driven by the observation of the quantum Hall effect in 
graphene even at room temperature [11 , the results on 
the graphene electronic structure modifications due to va- 
cancies in the presence of magnetic field are much sparser 

In this paper we consider the effects of vacancies in 
graphene lattices when a perpendicular magnetic field 
is applied, showing that vacancies give origin to extra 
(localized) states between Landau levels (LLs). When 
a sufficiently high density of vacancies is considered, we 
verify that these states can couple to each other, produc- 
ing additional bands with non-negligible density of states 
(DOS) between LLs. This problem has not been exhaus- 
tively explored in conventional semiconductor-based two- 
dimensional electron gases either. The line shape of Lan- 
dau bands, as well as the experimentally observed inter- 
LLs DOS [13], is a problem posed since the very begin- 
ning in the history of the quantum Hall effect, but many 
issues still remain without a conclusive answer [El |15] . 
The first available experimental observations of the LLs 
DOS in graphene [161 Hi] seem indeed to show a non- 



vanishing background DOS between LLs and, for the 
measurements on the surface of highly oriented pyrolytic 
graphite, unidentified DOS structures between LLs were 
reported [16 . 

The coupling between the wave functions of two va- 
cancies is verified to be determined by the distance be- 
tween them and by the magnetic length tuned by the 
magnetic field. When considering densities of vacan- 
cies, we introduce the vacancies in a periodic pattern, 
which is equally spaced in the graphene lattice. Strik- 
ingly, from the investigation of the localization proper- 
ties of the additional bands formed between the LLs due 
to these vacancy-patterns, we envisage that they could 
lead to interesting features in magneto-transport mea- 
surements, such as extra Hall resistance plateaus. A ran- 
dom distribution pattern is expected for vacancies seen 
as natural disorder sources or formed after irradiation 
with ion-beams. However, due to the advantageous elec- 
tronic properties of graphene, such as room-temperature 
ballistic transport [2, 3 , microstructuring efforts have ex- 
perienced a fast recent development after initial concept 
devices. Two groups have recently reported Coulomb 
blockade fingerprints in graphene quantum dots with typ- 
ical dimensions in the tens of nanometers range JHl |19] . 
Furthermore, the first attempts of tailoring graphene 
sheets with scanning probe microscope lithography have 
also been reported in the literature [20l|2T]. Therefore, 
envisaging new structures is a widespreading goal for 
graphene based applied investigations. Rethinking previ- 
ously tested structures and geometries on a different re- 
search context has always been a motivation for several 
groups. A recent example is the concept of an ant idot 
lattice on graphene [22 . Taking advantage of these char- 
acteristics, together with the fact that microstructured 
graphene may lead to a working device without capping 
layers, an even further step could be proposed: vacancies 
may be thought not as defects, but may be introduced on 
purpose by following a previous design. From this point 
of view, vacancies are the minimal possible inner edges 
or, in other words, antidot units. 



II. VACANCIES IN THE LATTICE MODEL 

The results shown here are obtained within a tight- 
binding model with nearest neighbors hoppings in a hon- 
eycomb lattice, a framework that adequately grasps the 
electronic properties of graphene monolayers. We cal- 
culate the energy spectrum of graphene layers with the 
mentioned vacancies in the presence of perpendicular 
magnetic fields. In this model, we use the following 
Hamiltonian: 



Furthermore, there are experimental evidences that the 
atomic positions of the carbon atoms in graphene do not 
reconstruct around a mono- vacancy [10]. In our model, 
the vacancies are introduced in graphene rectangular lat- 
tices of up to 8.000 atomic sites, which are repeated by 
means of periodic boundary conditions. 



III. VACANCY STATES BETWEEN LANDAU 
LEVELS 
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where q is the fermionic operator on site i. The mag- 
netic field B is included by means of a Peierls' sub- 
stitution in hopping parameter {t^2.7eV for graphene): 
(j)ij = 27r{e/h) p A-d\ . In the Landau gauge, (j)ij =0 
along the x direction and (j)ij = ±7r(x/a)<l>/<l>o along the 
TV direction, with ^/^o = Ba^^/?>e/{2h) (a=2.46A is 
the lattice constant for graphene). On-site disorder is in- 
troduced for the sake of LL bands broadening in the DOS 
calculations [23] : white- noise fiuctuations are included by 
sorting uncorrelated orbital energies within £i < \W/2\. 

We start by dealing with graphene systems containing 
isolated monovacancies or divacancies, as depicted in Fig. 
1. Different colors for the carbon atoms represented in 
Fig. 1 indicate the two different sublattices. Vacancies 
are simulated by setting the hopping parameters to zero 
and the on-site energy at the defect site equals to a large 
value outside the energy range of the DOS. The atomic 
relaxation contributions to the atomic geometry of the 
mono and divacancies in graphene are not considered. 
Reconstruction and rebounding around vacancies do not 
modify qualitatively the results we showed in this work. 
The main ingredient here is that vacancies are deep de- 
fects showing states that are strongly localized around 
them, as expected from numerical simulations [8] [24] . 




FIG. 1: (Color online) Representation of two typical va- 
cancies in the graphene lattice: Single vacancy (also called 
mono-vacancy) and divacancy (corresponding to two neigh- 
bor atomic vacancies). 



The complete electronic spectrum as a function of mag- 
netic fiux (the well-known Hofstadter-type spectrum) ob- 
tained for the perfect honeycomb lattice (no disorder) is 
depicted in Fig. 2(a). For the purpose of the present 
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FIG. 2: (Color online) (a) Energy spectrum as a function of 
magnetic flux for the hexagonal lattice. Rectangular region 
identifies the continuum limit of the spectrum, correspond- 
ing to the LLs of a monolayer graphene. Bottom: Energy- 
magnetic flux spectra for (b) a perfect lattice, (c) lattice with 
a single- vacancy and (d) lattice with a divacancy. Continuum 
lines show the energy dependencies with y/W]^o (i.e., with 
\/B) of the n=0,=bl,±2 graphene LLs, for guidance. Vacan- 
cies in the lattice introduce states with energies between the 
LLs, which are shown as open circles, for both vacancy types 
considered. 



work, we focus on the low-flux limit, which is identifled 
by the rectangular area and showed in more detail in 
Fig. 2(b). In this limit, with energy window around the 
Dirac point, the LLs are well deflned. Superposed to the 
spectra of Figs. 2(b), (c) and (d), continuum lines indi- 
cate (for guidance) the energy dependencies with VnB 
of the graphene LLs with indices n=0,±l,±2. For higher 
magnetic fluxes, the lattice effects are already seen as a 
broadening of the highly degenerate LLs and a deviation 
of n 7^ LLs from the continuum lines. 

In Fig. 2(c) one can observe the evolution of the states 
introduced by the absence of an atom (single vacancy) in 
the graphene sheet: a defect state appears in each energy 
range between consecutive LLs (open circles). Lowering 
the magnetic field brings the defect states to the cen- 
ter of the spectrum (zero mode energy) as already dis- 
cussed in comprehensive studies in the absence of mag- 
netic fields [8] . We observe that the energies of the states 
introduced by the vacancy are increased with magnetic 
fiux, however, these two states do not follow an energy 
dependency with V^ (followed by the LLs) nor have a 
linear increase with B. This result can be interpreted 
as the evolution of an edge state. Indeed, a vacancy is 
a minimal inner zigzag edge. The energies of external 
edge states are lowered with magnetic field, while states 
from inner edges should follow the opposite fate [25 . A 
divacancy is observed to introduce two states in the en- 
ergy ranges between consecutive LLs [open circles in Fig. 
2(d)]. These two states, as for the previous case, fol- 
low a different energy dependence from the LLs. From 
the point of view of an edge picture, a divacancy mixes 
zigzag and armchair characters. The inclusion of relax- 
ations around the vacancies is not expected to modify 
qualitatively the physical picture shown here. 



IV. TWO CLOSE VACANCIES: COUPLING OF 
THE LOCALIZED STATES 

If two single vacancies are introduced far from each 
other in the lattice, two degenerated states are intro- 
duced in each inter-LL energy window at the same en- 
ergy positions shown in Fig. 2(c) for only one va- 
cancy. However, a different scenario shows up if the 
two vacancies are introduced close to each other as re- 
vealed in the energy-magnet ic-fiux spectra shown in Figs. 
3(b-f). Five inter- vacancies distances are considered: 
D=5acc, D=6acc, D=7acc [schematically represented in 
Fig. 3(a)], D=Saccj and D=9acc^ where acc=l-42A is the 
carbon-carbon distance, i.e., the atomic spacing between 
two neighbor carbons in the graphene lattice. 

Under the main features identified in the spectra of 
Fig. 3 are the very energy splitting of the defect states 
and the evolution of this splitting with magnetic field. 
One expects that these levels split if there is a hybridiza- 
tion between the states localized around the vacancies. 



For an inter- vacancy distance of D = 9acc [Fig. 3(f)], 
the splitting observed is becoming small and the two 
vacancy-related states almost collapse in energy at high 
magnetic fields. Hence, one can see in Fig. 3 that the 
energy splitting decreases with increasing vacancy dis- 
tances or increasing magnetic field. These main features 
observed are robust independently of whether the two 
vacancies are situated in the same sublattice or in differ- 
ent sublattices. In fact, for the inter- vacancies distances 
D=5acc^ D=7acc, and D=Sacc shown here, the two va- 
cancies are in different sublattices, while for the cases 
D=6acc and D=9acc^ the two vacancies are located in 
the same sublattice. The zero mode states (in the ab- 
sence of magnetic field) associated to single vacancies are 
localized around the absent atoms [8] . In the presence of 
a magnetic field, this localization could be seen from the 
point of view of a different scale, given by the magnetic 
length lB = \/h/eB ^ 2b7 k/ ^B{Tesla)] increasing the 
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FIG. 3: (Color online) (a) Representation of two close va- 
cancies in the graphene hexagonal lattice, with a separation 
between them of D = 7acc, where ace = 1.42A is the dis- 
tance between two neighbor carbon atoms. Bottom: Energy- 
magnetic-flux spectra for lattices with two close vacancies, 
separated by increased distances: (b) 5acc, (c) 6acc, (d) 7acc, 
(e) Sttcc, and (f) 9acc- Two states due to the two vacancies 
appear between LLs (open circles) and the energy splitting 
between these states increases for decreasing distance D or 
decreasing magnetic fields. 



magnetic field will then lead to a spatial shrinking of the 
wave function. This shrinking leads to a quenching of 
the hybridization of the vacancy states, which will start 
to be effective only for Ib ^ D/2. 



V. VACANCY SUPERLATTICES 



A. Extra bands in the DOS 



We then consider higher vacancy concentrations, in- 
troduced in a periodic pattern. As we discuss in Sec. I, 
having in mind a promising and bearable scenario of mi- 
croscopic manipulation, a vacancy lattice - the ordered 
counterpart of a disordered network - is worth to be con- 
sidered. In Fig. 4(a) the DOS of a graphene lattice with- 
out vacancies is compared with the DOS of a lattice con- 
taining a rectangular arrangement of 12 vacancies, cor- 
responding to a vacancy concentration of x=0.21%. The 
dimensions of the lattices considered here in Fig. 4 are 
96x60=5760 atomic sites, arranged in 96 zig-zag chains, 
each of them containing 60 sites. Here a slight white- 
noise on-site disorder is also introduced {W/t = 0.1) in 
order to introduce a broadening in the highly degenerate 
LLs, as well as to permit an unambiguous calculation of 
the localization properties. 

The focus of the present work is on the extra bands 
created between the LLs [clearly shown in Fig. 4(a)]. 
For the magnetic flux considered, $/<^o = 2/96 ^ 0.02, 
the bands due this density of vacancies are already com- 
parable to the bulk LLs. It should be kept in mind that, 
for a fixed concentration of vacancies, the relative weight 
of these extra bands further increases with diminishing 
magnetic fields, since the number of states in each extra 
band is equal to the number of vacancies in the lattice, 
while the LL degeneracy is proportional to the magnetic 
field. In fact, the LL degeneracy is given by Ntotai x ^/^o^ 
the total number of atomic sites considered times the 
magnetic flux. Another point to be noticed is that, for 
inter-vacancies distances corresponding to this vacancy 
concentration of x=0.21%, the associated energy levels 
are still degenerated (when the white-noise disorder is 
not included). This can be verified by inspecting that 
the widths of the LLs and of these extra-bands are basi- 
cally the same, induced by the white- noise disorder. The 
higher concentration showed here, of 0.42% [in Fig 4(b)], 
imposes a vacancy lattice with distances between vacan- 
cies of I^ ~ ISttcc in both lattice directions. Therefore, 
vacancy (antidot) lattice effects are expected in the elec- 
tronic spectrum for magnetic fluxes at which the "bulk" 
graphene LLs are still well defined. This indeed is the 
origin of some non-trivial modulation observed for the 
DOS between the higher LLs (not shown here). 
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FIG. 4: (Color online) (a) Density of states for a graphene 
lattice without vacancies compared with a lattice containing a 
vacancy concentration of x=0.21% shifted vertically for clar- 
ity. The formation of extra bands between LLs, due to the 
vacancies, is evident. LLs are broadened by a white- noise 
disorder with W/t=0.1 and $/$o = 2/96 ^ 0.02. (b) DOS 
and participation ratio (PR) for states within the zero-energy 
LL and the first vacancy-related band. Increasing vacancy 
concentration (x) leads to a modulation of the PR within the 
extra band similar to that of the LL. 



B. Delocalization within vacancy bands 

We inspect the localization properties of the states 
within the extra vacancy bands by means of the par- 
ticipation ratio (PR) [26] , 



PR = 



N 



^EtilV'^ 



(2) 



where ipi is the amplitude of the normalized wave func- 
tion on site i, and N is the total number of lattice sites. 
The PR gives the proportion of N sites over which the 
wave function is spread. In this way, the PR for a local- 
ized state will tend to zero in the thermodynamic limit. 
In Fig. 4(b), attention is focused only on the Landau 
bands n=0 and n=+l. We see from the PR calculation 



that if the concentration of vacancies is low {x = 0.21%), 
ah the states in the vacancy bands are locaUzed. On the 
other hand, if the intervacancy distances diminish, at cer- 
tain magnetic-flux range (roughly when Ib > D/2) these 
vacancy states are effectively coupled and eventually be- 
come delocalized. Strikingly, the PR for a vacancy con- 
centration of X = 0.42% shows a strong modulation for 
the vacancy related band, similar to the ones for the Lan- 
dau bands [23]. This behavior signalizes that a vacancy 
lattice could be identified in magnetotransport measure- 
ment, possibly by extra quantum Hall plateaus. 



VI. DISCUSSION AND CONCLUSIONS 

For both vacancy concentrations shown in Fig. 4, all 
the vacancies were introduced in the same sublattice. 
However, we emphasize that we also tested configura- 
tions (not shown here) where the vacancies are equally 
distributed between the two sublattices (zero degree of 
uncompensation), and exactly the same effect seen in Fig. 
4 is observed, i.e., a vacancy-band is formed in the DOS 
and, after a critical concentration, the PR is modulated, 
indicating a delocalization of the states in the middle 
of the vacancy-band. In the absence of magnetic fields, 
it was shown ^ that for lattices containing vacancies 
with an increasing degree of uncompensation between 
sublattices, an increasing energy gap appears in the DOS 
around zero energy, with localized states in the middle of 
the gap. However, the magnetic field introduces a highly 
degenerate Landau level exactly at zero energy, which 
substantially modifies the scenario. Having in mind that 
for the concentrations and magnetic fluxes considered no 
changes are noticed due to sublattice uncompensation, 
we have that here, in the same way that for the case of 
two vacancies, the effect we are describing does not show 
dependence on the degree of uncompensation between 
sublattices. Nevertheless, the uncompensation effects on 
the limit of B^O, i.e, on the transition from well-defined 
LLs to B=0, deserve further investigations. 

We would like to mention that comparisons with exper- 
imental conditions are not hindered by the drastic lim- 
itations imposed by the periodic boundary conditions, 
which are further constrained by the commensurability 
related to the magnetic flux. For the largest lattice 
used in the results shown here (with 8000 sites being 
400 zig-zag chains, each with 20 sites in the results from 
Figs. 2 and 3), the lowest possible magnetic flux is still 
^/^o = 2/400 = 0.005, corresponding to B ^ 400T, 
a value that is one order of magnitude above attain- 
able experimental setups [27]. Intervacancy coupling, the 
main underlying cause of the present findings, is signif- 
icant whenever Ib > D/2. This condition implies, for 
B = 40T, an intervacancy distance D < 80A, which 
corresponds to a vacancy concentration of the order of 
x=0.1%. An intervacancy distance of SOAis also much 



lower than the estimations of both: mean-free path and 
coherence length in graphene samples. 

It should be emphasized that our results are not qual- 
itatively dependent on the precise position in the energy 
of the vacancy-related levels. As stated in Sec. II, the 
basic ingredients for building "vacancy molecules" or "va- 
cancy superlattices" are wave functions localized around 
the absent atoms that can be tuned by the magnetic field. 
The only condition in these additional levels is that they 
lie far apart in energy from the LLs, a very robust condi- 
tion. Besides that, an important fingerprint is the qual- 
itatively different evolution of the states with magnetic 
field, deviating from the usual y/B dependence for lateral 
LLs. 

In summary, our main findings are listed below, (i) Va- 
cancies in graphene introduce states between LLs, which 
tend to higher energies with increasing magnetic fields, 
consistent with the behavior of inner edge states. These 
states could be assessed by the usual switching of the 
magnetic field or gate voltages in quantum Hall mea- 
surements, but are not expected to reveal any fingerprint 
in these measurements because of their strong localized 
character, (ii) However, sufficiently near vacancies will 
lead to the coupling of these localized states, leading to 
"vacancy molecules" . The bonding between these vacan- 
cies can be tuned by the magnetic field, (iii) A periodic 
vacancy lattice introduces well-defined extra bands be- 
tween the LLs with a non- negligible DOS. They do not 
affect the transport properties only at low concentrations. 
However, in the low magnetic field limit or higher concen- 
trations, when vacancy states can be effectively coupled, 
vacancy superlattice delocalized states will cross the sys- 
tem and directly associated magneto-transport features 
can be expected. 
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